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Quantum gravity effects of zeroth order in the Planck constant are investigated in the framework 
of the low-energy effective theory. A special emphasis is placed on establishing the correspondence 
between classical and quantum theories, for which purpose transformation properties of the fi°-order 
radiative contributions to the effective gravitational field under deformations of a reference frame 
are determined. Using the Batalin-Vilkovisky formalism it is shown that the one-loop contributions 
violate the principle of general covariance, in the sense that the quantities which are classically 
invariant under such deformations take generally different values in different reference frames. In 
particular, variation of the scalar curvature under transitions between different reference frames is 
calculated explicitly. Furthermore, the long-range properties of the two-point correlation function of 
the gravitational field are examined. Using the Schwinger-Keldysh formalism it is proved that this 
function is finite in the coincidence limit outside the region of particle localization. In this limit, 
the leading term in the long-range expansion of the correlation function is calculated explicitly, 
and the relative value of the root mean square fluctuation of the Newton potential is found to be 
1/V2. It shown also that in the case of a macroscopic gravitating body, the terms violating general 
covariance, and the field fluctuation are both suppressed by a factor 1/N, where N is the number 
of particles in the body. This leads naturally to a macroscopic formulation of the correspondence 
between classical and quantum theories of gravitation. As an application of the obtained results, 
the secular precession of a test particle orbit in the field of a black hole is determined. 
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In this Chapter, quantum properties of the gravitational fields produced by quantized matter will be considered 
in the framework of the low-energy effective theory. This means that our interest will be in quantum aspects of the 
gravitational interaction at energy scales much lower than the Planck energy 



From the practical point of view, this does not present a limitation, since all phenomena observed in the Universe so 
far fall perfectly into this category. In view of the extreme smallness of the Newton constant G, these phenomena are 
well described by the lowest-order Einstein theory. On the other hand, the results obtained within the lowest-order 
approximation are model-independent, i.e., they play the role of low-energy theorems. This universality allows one to 
draw a number of important conclusions concerning the synthesis of quantum theory and gravitation. 

An essential part of investigation of this synthesis consists in establishing a correspondence between classical and 
quantum theories. The formal rules of this correspondence are contained in the Bohr correspondence principle which 
gives a general recipe for the construction of operators for various physical quantities. Informally, the problem of 
correspondence is in elucidating actual conditions to be imposed on a system to allow its classical consideration. 
Identification of these quasi-classical conditions constitutes an integral part of interpretative basis of quantum theory. 
Right from this point of view the low-energy behavior in quantum gravity will be considered below. 

Two characteristics of the gravitational field produced by a quantized system will be of our main concern in 
this Chapter - the mean value of the gravitational field, and its correlation function. The results of independent 
investigation of these quantities will lead us to one and the same macroscopic formulation of the quasi-classical 
conditions in quantum gravity. 



I. INTRODUCTION 




*A chapter contributed to "Progress in General Relativity and Quantum Cosmology Research" (Nova Science Publishers). 
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Turning to a more detailed formulation of our approach, let us note, first of all, that the low-energy condition stated 
above implies a similar condition on the characteristic length scale L of the process under consideration, L 3> lp, 
where lp is the Planck length 



In other words, our interest is in the long-range properties of the gravitational field produced by a given system. As far 
as this system can be treated classically, lp is the length scale characterizing quantum properties of its gravitational 
interactions, because lp is the only parameter entering the theory of quantized gravitational field through the Einstein 
action 



In the presence of quantized matter, however, another parameter with dimension of length comes into play, namely 
the gravitational radius 



This parameter appears, of course, already in classical theory. The question of crucial importance is whether r g has 
an independent meaning in quantum domain, representing a scale of specifically quantum phenomena. This question 
may seem strange at first sight, as r g does not contain the Planck constant h, an inalienable attribute of quantum 
theory. However, well known is the fact that gravitational radiative corrections do contain pieces independent of h. 
In the framework of the effective theory, they appear as a power series in tv/V, just like post-Newtonian corrections 
in classical general relativity. This fact was first clearly stated by Iwasaki The reason for the appearance of H° 
terms through the loop contributions is that in the case of gravitational interaction, the mass and "kinetic" terms in 
a matter Lagrangian determine not only the properties of matter quanta propagation, but also their couplings. Thus 
the mass term of, e.g., scalar field Lagrangian generates the vertices proportional to 



containing inverse powers of h. Naively, one expects these be cancelled by ft's coming from the propagators when 
combining an amplitude. One should remember, however, that such a counting of powers of h in Feynman diagrams 
is a bad helpmate in the presence of massless particles. Virtual propagation of gravitons interacting with matter 
field quanta near their mass shells results in a root non-analyticity of the massive particle form factors at zero 
momentum transfer (p). For instance, the low-energy expansion of the one- loop diagram in Fig. |2{a) begins with 
terms proportional to (— p 2 ) -1 / 2 , rather than integer powers (or logarithms) of p 2 . It is this singularity that is 
responsible for the appearance of powers of r g /r in the long-range expansion of the loop contributions. The question 
we ask is of what nature, classical or quantum, these pieces are. This is precisely the question of correspondence in 
quantum gravity. 

Suppose that the matter producing gravitational field satisfies the usual quantum mechanical quasi-classical con- 
ditions, e.g., consider sufficiently heavy particles. Then the quasi-classical conditions for the gravitational field can 
be inferred from the requirement that the mean value of the spacetime metric, (g^), coincides with the classical 
solution of Einstein equations, corresponding to the same distribution of gravitating matter. Practically, the most 
convenient way of looking for these conditions is to compare transformation properties of the quantities involved under 
deformations of the reference frame, thus avoiding explicit calculation of the expectation values. The latter point of 
view takes advantage of the fact that the transformation law of classical solutions is known in advance. Hence, we 
have to check whether {g^ u ) transforms covariantly with respect to transitions between different reference frames. In 
other words, we have to consider the question of general covariance in quantum gravity. This is the way we follow in 
Sec. IIHI Alternatively, one can infer the quasi-classical conditions from the requirement of vanishing of field fluctua- 
tions. We take this rout in Sec. I1VI Some results of the Schwinger-Keldysh and Batalin-Vilkovisky formalisms, used 
in our investigation, are summarized in Sec.[H] Section Ivl contains concluding remarks. 

Condensed notations of DeWitt are in force throughout this chapter. Also, right and left derivatives with 
respect to the fields and the sources, respectively, are used. The dimensional regularization of all divergent quantities 
is assumed. 




(1) 




(2) 
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II. PRELIMINARIES 



Before going into detailed discussion of the question of general covariance in quantum gravity, let us describe the 
general setting we will be working in. 



A. Frame of reference and interacting fields. 

First of all, we should set a frame of reference, i.e., a system of idealized reference bodies with respect to which the 
4-position in spacetime can be fixed. Let us assume, for definiteness, that the frame of reference is realized by means 
of an appropriate distribution of electrically charged matter. For simplicity, the energy- momentum of matter, as well 
as of the electromagnetic field it produces, will be assumed sufficiently small so as not to alter the gravitational field 
under consideration. The 4-position in spacetime can be determined by exchanging electromagnetic signals with a 
number of charged matter species. The electric charge distributions a a of the latter are thus supposed to be in a 
one-to-one correspondence with the spacetime coordinates x^, 



where index a enumerates the species. The <J a {x) will be assumed smooth scalar functions. Physical properties of the 
reference frame are determined by the action S a which specific form is of no importance for us. 

Next, let us consider a system of interacting gravitational and matter fields. The latter are arbitrary species, bosons 
or fermions, self-interacting or not, denoted collectively by <f>i, i = 1, 2, k. Dynamical variables of the gravitational 
field are h^ v = g^ v — f]^. Dynamics of the system is described by the action S g + S<j» where S<f, is the matter action, 
and S g is given 1 by Eq. (0). 

The total action S = S g + S<p + S a is invariant under the gauge transformations 

$V' = Cdah^ + (rif^a + h^dvC + (r] va + h va )d^ a = G^ a , (3) 
S4H = GfU, (4) 
Sa a = cr Q:Q £ Q (5) 

The generators Gn U , Gi span the closed algebra 



liv a\ [iv ^erA J 7 M 

Gt k G{-G^Gt = f a %G], (6) 



where the "structure constants" f a ^ 7 are defined by 

f a0 ^aVp = Had^ - Va d a ^ . (7) 
Let the gauge-fixing action be written in the form 

S gf = (f o - \^Cpc}j n a , (8) 

where F a is a set of functions of the fields /i M „, fixing general invariance, 7r Q auxiliary fields introducing the gauge, 
and C Q/ g a non-degenerate symmetric matrix weighting the functions F a ; the particular choice £ Q| g = $,r] a p corresponds 
to the well-known Feynman weighting of the gauge conditions. Introducing the ghost fields c Q , c Q , and denoting all 
the fields collectively by <&, we write the Faddeev-Popov action Q 

S FP m =S + S gl + cf 3 Ff u G%c a . (9) 



1 Our notation is Rp V = -R a Ma „ = doY^ — ■ ■ ■, R = R^g^" ', g = detg^, g^v = sgn(+, — , — , — ), ?7 M „ = diag{+l, —1, —1, —1}. The 
Minkowski tensor r] is used to raise and lower tensor indices. The units in which h = c = 16ttG = 1 are chosen in what follows. 
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5fp [3>] is invariant under the following Becchi-Rouet-Stora-Tyutin (BRST) transformations £| 





slot _ \ 

— ^^iV C Oi^l 


0<Pi 


slot, _ \ 
— Kj i C a A , 


S<7a 


— Oa.aC A, 


<5c 7 


2 J 7 c 


5c a 


= 7T Q A , 


Sn a 


= o, 



( cgA , 

(10) 

where A is a constant anticommuting parameter. Finally, the generating functional of Green functions has the form 
Z[J,K] = J'D<S>exp{i(j:{<P,K}+p a c a +p a c a +t^h^+ffa+ S a a a )}, (11) 

where 

£[$, X] - S FP [$] + + 9 *Gf c* + rV Q , Q c Q - y / a/3 7 c QC/3 + n Q 7r a , 

{i, j, s, (3, (3, 0} = J ordinary sources, and {fc, g, r, I, n, 0} = K the BRST-transformation sources for the fields 
{h, <f>, a, c, c, 7r} = $, respectively, and integration is carried over all field configurations satisfying 

$ ± — > , for £ — > q=oo , 

where the superscripts + and — denote the positive- and negative-frequency parts of the fields, respectively. The 
if-sources are introduced into Z for the purpose of studying gauge dependence of observable quantities. 



B. Effective fields and correlation functions 



The central objects of our investigation in the sections below will be the expectation value of the gravitational field 
produced by an elementary particle, and its correlation function. They can be obtained from the one- and two-point 
Green functions of the gravitational field, respectively, 



SW[J, K] 



8W V 
S 2 W\J,K] 



J\j=0 

k=o 



J\j=0 
K = Q 



(12) 
(13) 



with the help of the well-known reduction formula applied to the external matter lines to go over onto the mass shell 
of the particle. In Eqs. I|12fl . (|13H . W is the generating functional of connected Green functions, 

W[J,K] = -i\nZ[J,K], 

and the symbol J\j means that the source j is excluded from the set J, 

However, under the Fcynman asymptotic conditions specified above, Eqs. (|12|) . I|13[) give the in-out matrix elements 
of operators, rather than the in-in expectation values we are interested in. As is well known, in order to find the latter, 
the usual Feynman rules for constructing the matrix elements must be modified. According to the so-called closed 
time path formalism of Schwinger and Keldysh 0,0] (modern reviews of this formalism can be found in Refs. 0,OJ), 
this amounts to using in Eqs. (|12fl . (|13|) the generating functional 

Z SK [J ± ,K ± ] = y"l)$_y"2)$ + exp{z(E[$ ± ,^ ± ] + p a c a + (3 a c a + iT V 

+ 3% + s a a a )}, (14) 



where 



n 



E[$±, K ± ] = S FP [<P + } - 5 FP [*_] + k^G«c a + q l G a lCa + r a a a , a c a - -f^c^ + n a n a 
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instead of 111(1 . Here the subscript + (— ) shows that the time argument of the integration variable runs from — oo to 
+00 (from +00 to —00). Integration is over all fields satisfying 

<I> + — > , for t — ► -00 , 

$ + = for t -> +00. (15) 

We do not distinguish the plus and minus field components in the source terms explicitly, implying that summation 
over repeated Greek indices includes summation over ± as well as spacetime integration. For instance, 



etc. 



With the help of the new generating functional, the in- in expectation value of, e.g., the product h^ v (x)h a /3(x') can 
be written as 



(in I h^(x)h af 3 (x 1 ) | in) = -Z, 



S 2 Z SK [J ± ,K ± ] 



SK 



(16) 



J\j = 
K = Q 



It is seen that (in|/z Att ' (x)h a ^ (x')|in) is given by the ordinary functional integral but with the number of fields doubled, 
and unusual boundary conditions specified above. Accordingly, diagrammatics generated upon expanding this integral 
in perturbation theory consists of the following elements. There are four types of pairings for each component of the 
set $, corresponding to the four different ways of placing two field operators on the two branches of the time path. 
They are conveniently combined into 2x2 matrices according to 2 



D++ £+- 



For instance, pairings matrices for the gravitational and scalar matter fields read 



V i{4>i(x)^k{y))o i(T4>i(x)<j> k (y)) 



£> ua (x,y)= ( ^^WjvK^o i{hap{y)Kv{x))a \ 

where the operation of time ordering T (T) arranges the factors so that the time arguments decrease (increase) from 
left to right, and (-)o denotes vacuum averaging. The "propagators" Q^ap, ^ik satisfy the following matrix equations 

<5^ a/3 SW.A = -e*£, f g(2 ' , C = ^a + « (17) 

<8 ik S kl = -t6\ , & k = i , (18) 

0(Pi0q>k 

where denotes the free field part of the gauge-fixed action after the gauge introducing fields 7r Q have been 
integrated out, and e, i are 2x2 matrices with respect to indices +, — : 

1 \ / 1 



1/' 1 I -1 



2 Gothic letters are used to distinguish quantities representing columns, matrices etc. with respect to indices +, — 
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As in the ordinary Feynman diagrammatics of the S-matrix theory, the propagators are contracted with the vertex 
factors generated by the interaction part of the action, 5 mt [$], with subsequent summation over (+, — ) in the vertices, 
each "— " vertex coming with an extra factor (—1). This can be represented as the matrix multiplication of 2) M „ Q /3, £>ik 
with suitable matrix vertices. For instance, the <j) 2 h part of the action generates the matrix vertex QJ^^ which in 
components has the form 



Vku {x,y,z)=e stu 



8 3 S 



8h llv {x)8^i{y)8(j)k{z) 



where indices s, t, u take the values +, — , and e s t u is defined by 
4> line is represented in this notation by a column 



satisfying 



-e = 1 and zero otherwise. An external 



«"»-(» 

For future references, we give explicit expressions for various pairings of a single real scalar field 



D++(x,y) 



D-+{x,y) = i 



,—ik{x — y) 



(2tt) 4 TO 2 - k 2 - iO 

d A k 



D—(x,y) 



d 4 fc 



,—ik{x—y) 



(2tt) 4 k 2 - m 2 - iO ' 



(27T) 



6(k°)8(k 2 - m 2 )e- ik ^ , D+-{x,y) = D_ + (y,x) . 



(19) 



(20) 



C. Generating functionals and Slavnov identities. 

The functional £[<&±,lsr±] can be written as |Io) 

5if>[$±,K±]~ 



8$, 



where the so-called reduced action E r and the gauge fermion ^> are defined by 

£ r [$±, K±] = S + k^G%c a + q l G?c a + r a a a . a c a - y/ Q/3 7 c ttC/3 + n a ir a 



*[*±,#±] = fc^ V + Q*4H + r a o a + l a c a + n a c a + c a \F a - -tt% c 

Let us further simplify notation abbreviating k^h^ + q l <fii + r a a a + l a c a + n a c a = K and similarly for other sums 
of products of fields and sources or derivatives with respect to <&, J, K. Omitting also the subscript ±, the variation 
of £ under infinitesimal variation of the gauge conditions takes the form 



8Y\$,K] 



SAV[<f>,K] 8X[<f>,K] 



8$ 



8K 



The corresponding variation of the generating functional 



8Z SK [J,K]=i J J D$ + exp{i(E + J$)} 



8AV($,K) 8E(<f>,K) 



8<P SK 

Integrating by parts and omitting 8 2 Y,/8K8& ~ <5' 4 '(0) in the latter equation gives 



8Z SK [J,K]=iJ^ J J D$+A4<exp{i(£ + J*)}. 



(21) 
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Since E is invariant under the BRST transformation (|10fl , a BRST change of integration variables in Eq. I|14|) gives 
the Slavnov identity for the generating functional 

4^ = 0, (22) 
which allows one to rewrite Eq. H21[> in terms of the generating functional of connected Green functions as 

5Wsk[J,K\ = J (23) 
where (X) denotes the functional averaging of X . 



III. TRANSFORMATION PROPERTIES OF EFFECTIVE GRAVITATIONAL FIELD 

In this section, we will investigate the structure of the /i°-order contributions to the effective gravitational field, 
concentrating mainly on their transformation properties under deformations of the reference frame. The appearance 
of such contributions through the radiative corrections has its roots in the field properties of quantized matter sources. 
To make the field-theoretic aspect of the problem clearer, it is convenient to get rid of the purely quantum mechanical 
issues related to the quantum particle kinematics assuming matter quanta sufficiently heavy to neglect the position- 
velocity indeterminacy following from the Heisenberg principle. This assumption justifies the use of the term "particle 
at rest," and allows one to introduce a fixed distance r between the particle and the point of observation. 

The mean gravitational field produced by a massive particle is a function of five dimensional parameters - the 
fundamental constants h, G, c, the particle's mass m, and the distance r. Of these only two independent dimensionlcss 
combinations can be constructed, which we choose to be \ — lp/r and x — lc/r, where Iq = h/mc is the Compton 
length of the particle. The above assumption means that the gravitational field is considered in the limit x — > 0. For 
a fixed particle's mass, small values of a imply large values of r, leading naturally to the long-range expansion of the 
mean gravitational field. Accordingly, we will use where appropriate the terminology and general ideas of effective 
field theories [3 [HQ 



A. General covariance at the tree level. 



From the point of view of the general formalism outlined in the preceding sections, the classical Einstein theory 
corresponds to the tree approximation of the full quantum theory. The tree contributions to the expectation values 
of field operators coincide with the corresponding classical fields, and the effective equations of motion 

1 tT = 0. 



expressing the translation invariance of the functional integral measure, go over into the classical Einstein equations. 
Some of the tree diagrams representing their solution are shown in Fig. ^ The results of the preceding section allow 
one, in particular, to reestablish the general covariance of these equations. 

As was mentioned in the Introduction, coordinate transformations are replaced in the quantum theory by the field 
transformations. In particular, transformations of the reference frame are represented by variations of the fields at, 
induced by appropriate variations of the gauge conditions. Namely, it follows from Eq. 1)23(1 at the tree level that a 
gauge variation AF Q induces the following variations of the metric and reference fields 

&9liv 9fiv,a^ "F 9^a^ t u "F 9vot^ t ft > (24) 

Sa a = a a , a E a , S Q = (c a AV) , (25) 

where A "J is the corresponding variation of the gauge fermion. 

The functions 9^1,^0 a undergo the same variations (|24H . (|25H under the spacetime diffcomorphism 

as" + fls", 

with 5x^ = — S M . Let us consider any quantity entering the Einstein equations, for instance, the scalar curvature 
R. Under the above change of gauge conditions, the tree value of R measured at a point er° of the reference frame 
remains unchanged, 

MbWO)] - ^ + ^ = - ga- _ 0. (26) 
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Analogously, the tree contribution to any tensor quantity O^. . . (or O plA ••), for instance, the metric <7 M „ itself, calculated 
at a fixed reference point <t°, transforms covariantly (contravariantly) , as prescribed by the position of the tensor indices 
of the corresponding operator. This is the manifestation of general covariance of the classical Einstein theory in terms 
of quantum field theory. 



B. General covariance at the one-loop order. 

Let us now consider the transformation properties of the one-loop contributions. As in Sec. IIII Al we have to 
determine the effect of an arbitrary gauge variation on the value of the effective metric, and also on the functions cr a , 
i.e., on the structure of the reference frame. After that, the transformation law of observables, defined generally as 
the diffeomorphism-invariant functions of the metric and reference fields, can be determined in the way followed in 
Sec. IIII Al For definiteness, we will deal below with the scalar curvature R. Since we are interested in the onedoop 
contribution to the first post-Newtonian correction, we can linearize R in h^u '■ 

R = dfVh^ -Dh, h= rTh^ ■ (27) 

The transformation properties of R under variations of the weighting matrix £ aj g are considered in Sec. IIII B II and 
under variations of the gauge functions F a themselves in Sec. IIII B 21 



1. Dependence of effective metric on weighting parameters. 

Dependence of the effective fields on the weighting matrix (^ a p can be determined in a quite general way without 
specifying either the gauge functions F a , or the properties of gravitating matter fields. We begin with the classical 
theory in Sec. IIII B 1 a| and then consider the onedoop order in Sec. IIII B 1 bl 

a. The tree level. As we saw in Sec. IIII Al arbitrary gauge variations lead to the transformations of the classical 
fields, equivalent to the spacetime diffeomorphisms, and thus do not affect the values of R. In the particular case of 
variations of the weighting matrix, however, not only R, but also the effective fields themselves remain unchanged. 
This means that the structure of reference frames in classical theory is determined by the functions F a only. As this 
differs in the full quantum theory, a somewhat more detailed discussion of this issue will be given in this Section. 

To demonstrate (^-independence of the classical metric, let us first integrate the auxiliary fields 7r Q out of the 
gauge-fixing action JSJ), 

Sgf - 4 = iF^Fp , = 62 . (28) 

The classical equations of motion thus become 

s(s s + sL ) 



-T"", (29) 



where T^ v is the energy-momentum tensor of matter. Using invariance of the action S under the gauge transformations 

Sg^ = Cd a g^u + g^dui 01 + gvad^C = V M £„ , (30) 
and taking into account the "conservation law" V AI T A " / = 0, one has from Eq. (|29(l 

F a ^—^V*S(x-y) = 0. (31) 

The matrix M^{x,y) — §Fp/8g llv {x)W^ l S{x — y) is non-degenerate; its determinant A = det M% (x, y) is just the 
Faddeev-Popov determinant, and therefore A / 0. Hence, one has from Eq. £2) F a £, af) = 0, and, in view of non- 
degeneracy of £ Q/3 , F a — 0. The classical metric is thus independent of the choice of the matrix , and in particular, 
of the replacements F a — > A^Fp. One can put this in another way by saying that the weighting matrix has no 
geometrical meaning in classical theory. 

This differs, however, in quantum domain. The classical equations (|29|l are replaced in quantum theory by the 
effective equations 
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where T, g c ® vl and are the effective action, metric, and energy- momentum tensor of matter, respectively. In general, 
the fields g^, do not satisfy the gauge conditions F a = 0, and moreover, depend on the choice of the weighting matrix 
£a/3. ^"/'-independence is inherited only by the tree contribution. 

Dependence on the choice of the weighting matrix generally represents an excess of the gauge arbitrariness over the 
arbitrariness in the choice of reference frame; it is therefore a potential source of ambiguity in the values of observables. 
This dependence causes no gauge ambiguity of observables only if it reduces to the symmetry transformations. In 
other words, under (infinitesimal) variations of the matrix £ Q ^, the fields g^ must transform as in Eq. (|3(J|) 

with some functions S Q . It will be shown in the following Section that this is the case indeed. 

b. The one-loop level. Let us now turn to the examination of the gauge dependence of h° loop contribution to 
the effective gravitational field. This contribution comes from diagrams in which virtual propagation of matter fields 
is near their mass shells, and is represented by terms containing the root singularity with respect to the momentum 
transfer between gravitational and matter fields. In the first post-Newtonian approximation, the only diagram we 
need to consider is the one-loop diagram depicted in Fig. Ufa). As a simple analysis shows, other one- loop diagrams 
do not contain the root singularities, while many- loop diagrams are of higher orders in the Newton constant. For 
instance, the loop in the diagram of Fig. Hfb) does not contain massive particle propagators, and therefore, expands 
in integer powers (or logarithms) of the ratio a = —p 2 /m 2 . 

Let us show, first of all, that the Schwinger-Keldysh rules applied to the diagram in Fig.|2Ia) reduce to the ordinary 
Feynman rules of the standard S-matrix theory. According to the former, there are eight diagrams corresponding to 
the eight different ways of placing the three internal vertices of this diagram into the plus and minus branches of the 
time path (see Fig- - But all diagrams in Fig. [21 except the first two diagrams [2(a) and 2(b)] are zeros identically 
because of the energy-momentum conservation in the vertices, and the mass shell condition for external matter lines 
(massive particle cannot emit a real massless graviton). As to the diagram 2(b), its h° order contribution is also zero. 
Indeed, the time component of the momentum transfer 

p° = Vm 2 + (p + q) 2 - y/m 2 + q 2 = O (^) . 

To the leading order, therefore, p° is to be set zero, and the product of the two internal graviton propagators vanishes 
because the arguments of the corresponding ^-functions are opposite: 

D_ + (k)D + _(k+p) ~0(fc°)0(-fc°) = 0. 

The remaining diagram 2(a) contains only Feynman (causal) propagators, and therefore so do all other diagrams 
arising in the derivations of this section. Accordingly, we will omit the subscripts SK, ±, and denote the -D++ 
functions simply by D. 

It follows from Eq. 12: il) that under a variation A^f of the gauge fermion, variation of the effective gravitational field 

■ off _ sw 



has the form 

XK - , _ . _ , , , 

K = 

where J \ j means that the source j is excluded from J. We are interested presently in variations of the weighting 
matrix , therefore, 

or, integrating 7r" out, 

M>($,K) = ^C a pM fh FT (34) 

According to general rules, in order to find the contribution of a diagram with n external </>-lines, one has to take 
the nth derivative of the right hand side of Eq. (KS.''>[1 with respect to j z , multiply the result by the product of n 
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factors ei(q 2 — to 2 ), where q, are the 4-momentum and polarization of the external ^,-neld quanta, and set q 2 — to 2 
afterwards. The second term on the right of Eq. (|33l) is proportional to the source j % contracted with the vertex 
Gfc a . This term represents contribution of the graviton propagators ending on the external matter lines. Multiplied 



by (q — to ), it gives rise to a non-zero value as q 



only if the corresponding diagram is one-particle-reducible 



with respect to the 4>-\ine, in which case it describes the variation of hf£, under the gauge variation of external matter 
lines. It is well-known, however, that (^-operators must be renormalized 3 so as to cancel all the radiative corrections 
to the external lines . Therefore, this term can be omitted, and Eq. (|33[l rewritten finally as 



Sh^ - -C aP AZ 



7/ 



j\j=o 

K = Q 



(35) 



The one- loop diagrams representing the right hand side of Eq. I|35|l , which give rise to the root singularity, are shown 
in Fig. 0] Let us consider the diagram of Fig. Ufa) first. It turns out that this diagram is actually free of the root 

singularity despite the presence of the internal </)-line. The rightmost vertex in this diagram is generated by c a Fj , 

where F^ denotes the linear part of F 7 . The graviton propagator connecting this vertex to the <f>-lme can be expressed 
through the ghost propagator with the help of the equation 



(1),<tA 





D 



(7(0)/3 

f_LV 



I 



(36) 



where Dg is the ghost propagator defined by 



(37) 



respectively. Equation (|36[) is the Slavnov identity (122[) at the tree level, differentiated twice with respect to t^", (3 a . 
Using this identity in the diagram Fig. Ufa) we see that the ghost propagator is attached to the matter line through 

the generator D$ a . On the other hand, the action is invariant under the gauge transformations © - JSJ, 



S<j>i 1 5hav ^ 



(38) 



Differentiating Eq. (|38|l twice with respect to (j>, and setting h^ v = yields 



5S^ 5G1 SSi 2) SGI 



S 3 S d 



5<fii8(j)k $4>i 5<f>i5(pi 54>k 5h a \54ii8(j)k 
Taking into account also the mass shell condition 

45-4 = o, 



)%o<Pk 



G^ = 0. 



(39) 



h=0 



we see that under contraction of the vertex factor with the external and internal matter lines, the (/>-particle propagator, 
Dik, satisfying 



2o(2) 



S 2 S 



-D u = -5, 



i • 



is cancelled 



D 



S 2 S { I ] 



ki 



G«°i a = G^ 



(40) 



h=0 



3 One might think that the gauge dependence of the renormalization constants could spoil the above derivation of Eq. 1231 . In fact, this 
equation holds true for renormalized as well as unrenormalized quantities 1 1 1 H . 

4 The above discussion is nothing but the well-known reasoning underlying the proof of gauge-independence of the S-matrix 111! . 
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We conclude that the H° contribution of the diagram Fig. Efa) is zero. As to the rest of diagrams, they are all 
proportional to the generator D$ a . Thus, the right hand side of Eq. I|35|) can be written 

Shf u = G^ a E a + 0(h) , E a = ±Di( 07 Ae 5 (Fs) . 

Since E a are of the order G 2 , one can also write, within the accuracy of the first post-Newtonian approximation, 

5hf v = G%E a , (41) 

where G"^ are defined by Eq. © with /i M „ — > hfjH ■ 

We thus see that under variations of the weighting matrix, the effective metric does transform according to Eq. iKSl't . 
To determine the effect of these variations on the values of observables, one has to find also the induced transformation 
of the reference frame, i.e., of the functions a a . Obviously, the gauge variation of cr a 's is represented by the same set 
of diagrams pictured in Fig.0Jb,c,d), 5 with the only difference that the leftmost vertex (fiv) in these diagrams is now 
generated by a a , a c a instead of D^ u c a . Thus, under variations of the weighting matrix, the functions a a (x) transform 
according to 

6a a = a a<a E a , (42) 

where S"'s are the same as in Eq. 11411 . 

Equations (|41l) and 1|42|) are of the same form as Eqs. (|24[1 and (|25p. respectively, which implies that the value of 
any observable O is invariant under variations of the weighting matrix, 

SO[h eS (x(a°))\ = #Utf£ + = |^S Q - |° H" = . (43) 

Stiffi p dx^ dx a dx» v ' 

In particular, 

6R[g eS (x(a ))] = 0. 

Furthermore, any tensor quantity O a/ g... (or O m " ' ), calculated at a fixed reference point er°, transforms covariantly 
(contravariantly) , as prescribed by the position of the tensor indices of the corresponding operator. This is in accord 
with the principle of general covariance. 



2. Dependence of effective metric on the form of F a 

Having established the general law of the effective metric transformation under variations of the weighting matrix, 
let us turn to investigation of the variations of the functions F a themselves. 

According to the general equation , a variation AF a induces the following variation in the effective metric 



6h« = 6{daAFa) 



5k» v 



(44) 

J\j=o 

K = 



The general structure of diagrams representing the one-loop contribution to the right hand side of this equation is the 
same as before and given by Fig.0] Contribution of diagrams (b), (c), and (d) is again a spacetime diffeomorphism. In 
the present case, however, diagram (a) gives rise to a non-zero contribution already in the order ft P. Namely, it is not 

difficult to show that the combination AF^'^ D^ U<J \ cannot be brought to the form proportional to the generator 

G(°) . Note, first of all, that the variation of AF a v D^ux with respect to £ Q/3 is proportional to G^; in the highly 
condensed DeWitt's notation, 

5(AF^ 1} D) = AF^DiF^SZF^^D = AF^DF^S^DG^, 

where the Slavnov identity (|36|l has been used. Hence, without changing the H° part of diagram (a), Cap can be set 
zero, in which case Eq. {gjj gives f[ 1] D = 0. Suppose that AFP'"" D^x = X^G^ , or shorter, AF^D = XG^°\ 



5 The diagram of Fig. U a ) does not contribute in this case. 
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with some X. Then one has = AF^DF^ = XF^G^ = XM(h = 0). Since the Faddeev-Popov determinant 
dct M 7^ 0, it follows that X = 0. Thus, the argument used in the preceding section does not work, and the question 
is whether contribution of the diagram (a) can be actually represented in the form G"„S Q . 

The answer to this question is negative, as an explicit calculation shows. This will be demonstrated below in the 
simplest case of a single scalar field described by the action 

s<p = \ [d^xV^gig^d^d^ ~ m 2 2 } , (45) 



and linear gauge conditions 

F.. = n^ v B..h.... - ( 

where g is an arbitrary parameter. According to Eq. (|44|l . the ^-derivative of the effective metric is given by 



F 1 = rrdph^ - ( ^—J )d 1 h, h = rTh^ , ( a p = , (46) 



dhf v _ S(c a dF a /dg) 



dg Sk^ 



J\j=o 

K = 



(47) 



There are two diagrams with the structure of Fig.^a), in which the scalar particle propagates in opposite directions. 
They are represented in Fig. [SJ In fact, it is sufficient to evaluate cither of them. Indeed, these diagrams have the 
following tensor structure 

aiq^qu + a%{p il q v + p v q^,) + azP^Pv + (Wlnv, 

where a{, i = 1, 4, are some functions of p 2 . When transformed to the coordinate space, the second and third terms 
become spacetime gradients, hence, they can be written in the form G^H,,. As was discussed in the preceding 
sections, the terms of this type respect general covariance, therefore, we can restrict ourselves to the calculation of 
ai and 0,4 only. On the other hand, diagrams of Fig. 03 go over one into another under the substitution q — > — q — p 
which leaves 0,1,0,4 unchanged. Thus, we have 

d 3 q d 3 p a(q)a*(q + p) 
(2^)3(2^)3 ^/2e q 2e q+p 

= I^ ) (p,q) + I 3 ^ ) (P,Q), 

I 3 »l a) (P,q) = -^J^{lw a ^ 5 (q 1 +p 7 )(q S + ks)-m 2 ^ 

xD(q + k) -W praX q a {q X + k x ) - m 2 ^- D pT7TU] (k)- 



= / / „J 3 (o \s L Z e*%„(p, q) , Po = e q+p -e qi (48) 



2 2 J dg 

,-(k) {-(k ( - PC )S$ + Sftkv + 6*X} D xea p{k - P ) , 

~C\p,q) = ili a) {p,-q-p), 

where a(q) is the momentum space amplitude for the scalar particle, normalized by 



(2») 

/i is an arbitrary mass scale, e q = \J q 2 + m 2 , and e = 4 — d, d being the dimensionality of spacetime. Explicit 
expressions for the propagators 

D^a\(k) = — _ QiV^ukakx + Tloxkpkv)-^ 

+ (yiixokvkx + f)p\k v k a + T\ va kyk\ + Tj^xk^k^)-^ 

+ (3g 2 - 4g)k^k u k a k x ^ , (49) 
D a (k) - ^ g kak0 

- k 2 ~ 2 fc 4 > 
1 
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Calculation of l|48[l can be further simplified using the relation 

dFl D- F dD 
dg 1 dg ' 

which follows from FiD = 0, and noting that all gradient terms in the graviton propagators, contracted with the 4> 2 h 
vertices, can be omitted (see Sec. IIII B"TJ) . i.e., only the first line in Eq. I|49fl actually contributes. 

Upon extraction of the leading contribution, Eq. I|48|) considerably simplifies. Note, first of all, that to the leading 
order in the long-range expansion, one has £ q + P ~ £ q ~ to, and hence, po i=a 0. Next, take into account that the 
function a(q) is generally of the form 

a(q) = b(q)e- i ^, 
where x is the mean particle position, and b{q) is such that 



J d 3 q b{ q yi x = 



(50) 



for x outside of some finite region W around x = 0. In the long-range limit, b(q + p) may be substituted by 6(g) : 
This implies that we disregard spatial spreading of the wave packet, neglecting the multipole moments of the particle 
mass distribution. Hence, Eq. (|48|l can be written as 



d 3 q d?p 
2m J J (2tt) 3 (2tt) 3 



\b(q)\ 2 e-^ x -^I^(p,q), p° = 0. 



Let the equality of two functions up to a spacetime diffeomorphism be denoted by 
multiplications in Eq. I|48[). and omitting terms proportional to p M , one obtains 



(51) 

Then, performing tensor 



'k 1 



1 



1 



/3(a) ( D a ) ^ _j„< 
V» \PW V J ( 27r )4 fc 4 ( fc + p )2 m 2 _ ( fc + g) 2 

{Vv [^k i + 2PQ)(P-Q)(Q-m 2 ) 

+ (|p 2 (fc 2 + 2Q) + (k + p) 2 ™ 2 ) (k 2 + P)(Q - to 2 ) 
+k^k v [g(P 2 + k 2 m 2 )(Q - to 2 ) + 4£>PQto 2 - 3gPQ 
+ |p 2 (P - 2Q)(Q - to 2 ) + 2gQ 2 {Q - to 2 ) - gPrn 4 

+2{k+p) 2 (Q{Q - 2m 2 ) - P(Q - to 2 ) + to 4 )] 
+2{k„q v + Kq^){k+pf{Q - P){Q - to 2 ) 



-2q^q v {k+pf{k 2 + P)(Q-m 2 )} 
Introducing the Schwinger parameterization of denominators 

1 



1 f°° 

-^2=-J dycxp{yk 2 }, 



1 



k 2 + 2{kq) 
one evaluates the loop integrals using 



(k+p) 2 
dzexp{z[k 2 + 2(kq)]} 



Q = (kq), P=(kp). 

- I dxexp{x(k + p) 2 } . 
Jo 



(52) 



J d d kexp{k 2 (x + y + z) + 2¥ l (x PfJi + zq^ L )+p 2 x} 



x + y + z 



d/2 



exp 



9 9 9 

p z xy — m^z* 
x + y + z 



J d d k k a exp{fc 2 (a; + y + z) + 2W l (xp ll + zq^) +p 2 x} = 



x + y + z 



d/2 



exp 



2 9 9 

p xy — to z 
x + y + z 



xp 



zq Q 



x + y + z 
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etc., up to six fc-factors in the integrand. This calculation can be automated to a considerable extent with the help 
of the tensor package f° r the REDUCE system. Changing the integration variables (x, y, z) to (t, u, v) via 



t{l + t + u)v 



m 2 (l + atu) ' ^ rn 2 (l + atu) ' 
integrating v out, subtracting the ultraviolet divergence 6 

1 ffi 



u(l + t + u)v 



(l + t + u)v 
m 2 (l + atu) 



I; 



:i o)dl >^)=-T^V(-)' 

lQn 2 e \mJ 



<? M ^ + V^{P - 2m ) 24 



setting e = 0, omitting gradient terms, and retaining only the ^-contribution, we obtain 



///- 

OO /> OO 



+ 



16?r 2 
H 3 m 2 



dudt 



o Jo 
2# 2 



?7^g ( 1 1 



D 2 



1 - 



H s a \D 2 2D / 

1 fig lip + 4 



D = 1 + crai . 



1 

Dj ' a \D 3 D 1 
H = 1 + u + 1 . 



D 



(53) 



The root singularity in the right hand side of Eq. (|53|) can be extracted using Eqs. i|A3[l derived in the Appendix. 
Denoting 



oo />oo 



JO 



dudt H- n D~ m = J„ 



one has 



J- 



12 



4^ 



rroot 

Ji q — 



7 J 13 



Troot 
J 31 



7T' 

16 



/~~ Troot 

V " , J32 



3tt 2 

16Va ' 
_3tt 2 
32 



15tt 2 
128 



Substituting these into Eq. gives 



128Va 

Finally, restoring ordinary units, substituting Eq. I|54[l into Eq. I|51|l. and using the formula 

d 3 p e l P x 1 



(54) 



(2tt)3 \p\ - 2vr 2 r 2 ' f ~ ' 
we obtain the following expression for the ^-derivative of the G 2 S°-order contribution to the effective metric 01 



(55) 



( X ) ~ o„2„2 



dg dg dg 



Og 



6° 5° 

—o-(Q+ 1) - 



(56) 



A technicality must be mentioned here. By itself, the diagram of Fig. £2a) is free of infrared divergences. As a result of the BRST- 
operating with this diagram, however, some fictitious infrared divergences are brought into individual diagrams representing the right 
hand side of Eq. 1331 . This is because the vertex D"C a contains the term C a d a h^ v in which the spacetime derivatives act on the 
gravitational, rather than the ghost field. These divergences occur as u, t — + 00. They are proportional to integer powers of p 2 , and 
therefore do not interfere with the part containing the root singularity. Since these divergences must eventually cancel in the total sum 
in Eq. 1331 . they will be simply omitted in what follows. 
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The right hand side of this equation cannot be represented in the form (|32() . This result can be made more expressive 
by calculating the p-variation of the scalar curvature R in a given point er° of the reference frame 

A^cff, , SR ( dtf™ dh%»\ OR ^ 

an /,loop op 

og oq c 4 r 4 

7 

or 

W))] = gV 

OQ c 4 r 4 

Equations l|56|l . (|58|l express violation of general covariance by the loop corrections. 

Thus, despite their independence of the Planck constant, the post- Newtonian loop contributions turn out to be of a 
purely quantum nature. 

We are now in a position to ask for conditions to be imposed on a system in order to allow classical consideration 
of its gravitational field. Such a condition providing vanishing of the h° loop contributions can easily be found out by 
examining their dependence on the number of particles in the system. Let us consider a body with mass M, consisting 
of a large number N = M/m of elementary particles with mass to. Then it is readily seen that the n-loop contribution 
to the effective gravitational field of the body turns out to be suppressed by a factor 1/N n in comparison with the 
tree contribution. For instance, at the first post-Newtonian order, the tree diagram in Fig. ^b) is bilinear in the 
energy-momentum tensor (T^ v ) of the particles, and therefore proportional to (to • N) ■ (to • N) = M 2 . On the other 
hand, the post-Newtonian contribution of the diagram in Fig. G{a) is proportional to to 2 • N = M 2 /N, since it has 
only two external matter lines. 

Thus, we are led to the following macroscopic formulation of the correspondence principle in quantum gravity: the 
effective gravitational field produced by a macroscopic body of mass M consisting of N particles turns into corresponding 
classical solution of the Einstein equations in the limit N — > oo 19] . In particular, the principle of general covariance 
is to be considered as approximate, valid only for the description of macroscopic phenomena. 

The fi°-order loop contributions are normally highly suppressed. For the solar gravitational field, for instance, 
their relative value is m pro ton/-ZW© ~ 10~ 57 . However, they are the larger the more gravitating body resembles an 
elementary particle, and can become noticeable for a sufficiently massive compact body. 



C. Effective gravitational field of a heavy particle 

As an application of the obtained results, we will calculate the effective gravitational field of a particle with mass 
M in the first post-Newtonian approximation. 

The complete expression of the order h° for the spacetime metric is the sum of two pieces. The first is the tree 
contribution corresponding to the classical Schwarzschild solution 

ds 2 cl = (l - -) c 2 dt 2 - - r 2 (d6 2 + sin 2 9 dip 2 ), (59) 

^ r 

where 9, ip are the standard spherical angles, r is the radial coordinate, and r g = 2GM/c 2 is the gravitational radius 
of a spherically-symmetric distribution of mass M. The second is the one-loop post-Newtonian correction contained 



7 Another way to obtain this result is to introduce the sources tR and kR >a c a for the scalar curvature and its BRST- variation, respectively, 
into the generating functional 1111 . instead of the corresponding sources for the metric. Then Eq. 1471 is replaced by 

dR BS = 5(c a dF a /d e ) 
dp Sk j\j=o ' 1 ' 

if = o 

At the second order in G, the nontrivial contribution comes again from the diagram of Fig. UTal in which the lower left vertex is now 
generated by i?, iQ c Q . Thus, only the linear part of R gives rise to a non-zero contribution to the right hand side of Eq. 1571 . In other 
words, 8R[h eB ] '= SR eB , though generally R[h cB ] ^ R cS . 
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in the diagram (a) of Fig. Restoring the ordinary units, and using expressions (55), (65) of Ref. [3 for the vertex 
formfactors, and (|49|l with g = for the graviton propagator, 8 we obtain 

h^(p) = -^hf ( 3M V + ^ + 7M 2 ^) . (60) 
cry—p \ c P / 

Written down in the coordinate space with the help of the formulas l|55|l and 



d 3 p PtPk „ 4px _ 1 ( , 2x,x k 

Oik o 



(2tt) 3 |p| 3 27r 2 r 2 



equation l|6U|) gives, in the static case, 



Before adding the two contributions, however, one has to transform Eq. (|59|l which form is fixed by the requirements 
gn = 0, i = r, 9, tp; gee = — r 2 , to the DeWitt gauge 

F 1 = rfd^h^ - -d^h , Ca/3 = Vafi , (62) 

under which Eq. I|6U|I was derived. According to the general theorems about (^-independence of the fi°-order 
contributions, proved in Sec. IIII B II we can transform our expressions to the singular case <^ a p = instead of 
Cq/3 — VaP- Then the effective gravitational field will satisfy 

r ) ^d„hf 7 -^d 1 h efi = 0. (63) 

The t, 8, ^-components of Eq. (|63|l are already satisfied by the classical solution l|59|l . To meet the remaining condition, 
let us substitute r — ► /(r), where / is a function of r only. Then the t, 8, (^-components of Eq. (|46|l are still satisfied, 
while its r-component gives the following equation for the function f(r) : 

1 d ( r 2 f' 2 \ 2/ 2 1 d ( r„ f' 2 2f 2 \ n 



r 2 dr \1 — Tg/ f J r 2 2 dr \ f 1 — r g /f r 
where /' = df(r)/dr. 

Since we are interested only in the long-distance corrections to the Newton law, one may expand f(r)/r in powers 
of Tgjr keeping only the first few terms: 



f(r) 



l + d^ + c/ 7 " 9 ^ 



r \ r 

Substituting this into Eq. (|64J) . one obtains successively c\ — 1/2, C2 = 1/2, etc. Therefore, up to terms of the order 



r 2 /r 2 , the Schwarzschild solution takes the following form 



_ r 2[ 1 + ^ + ^|](^2 + sin 2 e ^ 2) ( g 5) 



Rewriting Eq. I|61l) in spherical coordinates and adding the tree contribution, we finally obtain the following expression 
for the interval 201 



ds 2 cS ee g^dx'dx" = (1 - c 2 dt 2 ( 1 + ?f - 11| 1 dr 2 



7r 



2 



r y \ r 



1 + t5 + tJ ( d6> + sin 61 ) • ( 66 ) 



r 4r 



8 Note the notation differences between Ref. Il5l and this Chapter. 
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In connection with this result it should be noted the following. Taking into account higher-order radiative post- 
Newtonian corrections will result in a further modification of the Schwarzschild solution. Since quantum contributions 
are of the same order of magnitude as those given by general relativity, this modification can lead to a significant shift 
of the horizon. In particular, the metric g c ^(r) may well turn out to be a regular function of r when all the ft loop 
corrections arc taken into account. 



IV. QUANTUM FLUCTUATIONS OF GRAVITATIONAL FIELD 

We turn now to another aspect of the long-range behavior in quantum gravity - quantum fluctuations of the 
gravitational fields produced by elementary systems. On various occasions, the issue of fluctuations has been the 
subject of a number of investigations (see Refs. 0, H2, |2^, 0, 113 j where references to early works can be found). 
It should be mentioned, however, that despite extensive literature in the area, only vacua contributions of quantized 
matter fields to the metric fluctuations have been studied in detail. At the same time, it is effects produced by real 
matter that are of special interest concerning the structure of elementary contributions to the field fluctuation. 

Before wc proceed to calculation of the correlation function, we shall examine its general properties in more detail. 
Namely, the structure of the long-range expansion of the correlation function, and the question of its gauge dependence 
will be considered. 



A. Properties of correlation function 

As defined by Eq. . the correlation function C^ap is a function of two spacetime arguments (x,x'). Of special 
interest is its "diagonal element" corresponding to coinciding arguments and Lorentz indices (x = x',[iv = a/3), 
which describes a dispersion of the spacetime metric around its mean value in a given spacetime point. However, it 
is well-known that this element is not generally well-defined because of the singular behavior of the product of field 
operators in the coincidence limit x — ► x' . This difficulty is naturally resolved when one takes into account the fact 
(realized long ago, see, e.g., Ref. j2(|) that in any field measurement in a given spacetime point, one deals actually 
with the field averaged over a small but finite spacetime domain surrounding this point. Thus, the physically sensible 
expression for the field operator is the following 



Jf^ = — / dt d 3 x h^(x,t). (67) 



1 

VT J T jy 

Respectively, the product of two fields in a given point is understood as the limit of 

2V<*/3 = J^dt J^dt 1 J^d 3 x J^d 3 x' h^(x,t)h a p{x',t') (68) 

when the size of the domain tends to zero. Finally, correlation function of the spacetime metric in this domain is 

G^a/3 = ("B^ap) - (ftju,)(lk a/ 9) = fyyp J J d^xd^X C^ al3 (x , x') . (69) 

Figure [5] depicts the tree diagrams contributing to the right hand side of Eq. I|16|) ■ The disconnected part shown in 
Fig. EUa) cancels in the expression for the correlation function, Eq. (|69J) . which is thus represented by the diagrams 
(b)-(h). 



1. Correlation function in the long-range limit 

As we mentioned in the beginning of Sec. IIIII the mean gravitational field produced by a massive particle is a 
function of five dimensional parameters - the fundamental constants ft, G, c, the particle's mass m, and the distance 
between the particle and the point of observation, r. Unlike the mean field, however, G^a/3 is a function of two 
spacetime arguments. On the other hand, we are interested ultimately in the coincidence limit of this function, which 
will be shown below to exist everywhere except the region of particle localization. In this limit, therefore, G^ va p 
depends on the same five parameters h,G,c,m,r. Assuming as before the particle sufficiently heavy [>c — ► 0), we 
conclude that the relevant information about field correlations is contained in the long-range asymptotic of C /JJ/Q ,^(r). 

To extract this information we note, first of all, that in the long-range limit, the value of G^ a p is independent 
of the choice of spacetime domain used in the definition of physical gravitational field operators, Eq. I|67|l . Indeed, 
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in any case the size of this domain must be small in comparison with the characteristic length at which the mean 
field changes significantly. In the case considered, this requires that F « r 3 . To the leading order of the long-range 
expansion, therefore, the quantity (h l/tu (x)h a ff(x')) appearing in the right hand side of Eq. can be considered 
constant within the domain. However, one cannot set x = x' in this expression directly. It is not difficult to see 
that the formal expression {h^ v {x)h a i3{x)) does not exist. Consider, for instance, the diagram shown in Fig. E^b). A 
typical term in the analytical expression of this diagram is proportional to the integral 

where is the 4- momentum of the particle, p^ the momentum transfer, and V^ va p a vertex factor combined of 's 
and momenta q,p, k. For a small but nonzero (x — x') this integral is effectively cut-off at large fc's by the oscillating 
exponent, but for x = x' it is divergent. This divergence arises from integration over large values of virtual graviton 
momenta, and therefore has nothing to do with the long-range behavior of the correlation function, because this 
behavior is determined by the low-energy properties of the theory. Evidently, the singularity of I^ V up{x — x\p) for 
x — > x' is not worse than ln(x^ — x' tl ) 2 /(x tl — x'^) 2 . Therefore, (m\h tiu (x)h a p(x')\m) is integrable, and G^ap given 
by Eqs. (EH is w eU defined. 

Our aim below will be to show that this singularity can be consistently isolated and removed from the expression 
for Ifivapix — x' ,p) and similar integrals for the rest of diagrams in Fig. EI without changing the long-range properties 
of the correlation function. After this removal, it is safe to set x = x' in the finite remainder, and to consider G^ a /3 as 
a function of the single variable - the distance r. An essential point of this procedure is that the singularity turns out 
to be local, and hence does not interfere with terms describing the long-range behavior, which guaranties unambiguity 
of the whole procedure. 

It should be emphasized that in contrast to what takes place in the scattering theory, the ultraviolet divergences 
appearing in the course of calculation of the in- in matrix elements in the coincidence limit are generally non-polynomial 
with respect to the momentum transfer. The reason for this is the different analytic structure of various elements in 
the matrix propagators £> Mi , ai a, 3D, which spoils the simple ultraviolet properties exhibited by the ordinary Feynman 
amplitudes. 9 Take the above integral as an example. Because of the delta function in the integrand, differentiation 
of Ifivapix — x' ,p) with respect to the momentum transfer does not remove the ultraviolet divergence of I^ U ap{^,p)- 
What makes it all the more interesting is the result obtained in Sec. II V 51 below, that the non-polynomial parts of 
divergent contributions eventually cancel each other, and the overall divergence turns out to be completely local. 

Next, let us establish general form of the leading term in the long-range expansion of the correlation function. In 
momentum representation, an expression of lowest order in the momentum transfer with suitable dimensionality is 
the following 

CUP) - G 2 m 2 ^| , (71) 



where 'J >JlVa p is a dimensionless positive definite tensor combined from r)^ and q^. 

Not all of diagrams in Fig. El contain contributions of this type. It is not difficult to identify those which do not. 
Consider, for instance, the diagram (h). It is proportional to the integral 

which does not involve the particle mass at all. Taking into account m 2 coming from the vertex factor, and (2e q )^ 1 ^ 2 ~ 



ypm from each external matter line, we see that the contribution of the diagram (h) is proportional to m, not m 2 . 
The same is true of all other diagrams without internal matter lines. As to diagrams involving such lines, it will be 
shown in Sec. HVBl bv direct calculation that they do contain contributions of the type Eq. Ij71|) . But prior to this 
the question of their dependence on the gauge will be considered. 



As is well known, the proof of locality of the S-matrix divergences relies substantially on the causality of the pole structure of Feynman 
propagators. This property allows Wick rotation of the energy contours, thus revealing the essentially Euclidean nature of the ultraviolet 
divergences. 
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2. Gauge independence of the leading contribution 



As in the case of mean gravitational field considered in Sec. Iffll the 0-component of the momentum transfer is to 
be set zero when calculating the leading term in the long-range expansion of the correlation function. This implies, 
in particular, that 6^ a ^ contains information about fluctuations in a quantity of direct physical meaning - the static 
potential energy of interacting particles. As such, it is expected to be independent of the choice of gauge conditions 
used to fix general covariance. More precisely, we have to verify that under variations of the gauge conditions, Q^ a p 
varies in a way that does not affect the values of observables built from it . 

As we saw in Sec. IIV A II the only diagrams contributing in the long-range limit are those containing internal 
massive particles lines. We will show presently that the gauge-dependent part of these diagrams can be reduced to 
the form without such lines. First, it follows form Eq. 1|17|) that the gauge variation of the graviton propagator satisfies 

«W = ^^Se^^x ■ (72) 

On the other hand, contracting Eq. (|17fl with G$ p yields 

Defining the matrix ghost propagator according to Eqs. (|17fl . (|37fl by 



one finds 



Taking into account that 



ffv* = CnGT®} > ( 73 ) 



and using Eq. {73J, we get 

*2W = (^""aW i fy) G ( °l 7 + G<f (pi i SF^V^ax) ■ (74) 

Let us assume for definitcness that the pair of indices [iv refers to the point of observation, and consider the first 
term in Eq. I|74|) . This part of the gauge variation of the graviton propagator is attached to the matter line through 

the generator G^j 1 . Suppressing all Lorentz and matrix indices except those referring to the <p 2 h- vertex, it can be 
written as 

The matrix vertex Q3 crA < I ' c is obtained multiplying S 3 S < f,/d(pi6(/)kSh„x by the matrix e s t u - Equation l|39[) then shows 
that the combination b (T A2J <TA:lfe can be written as 



where 



9l0(pk OCpi 0(piC 



— a 7 _ 



Finally, contracting b a \SO aX ' lk with the matrix J)fe m and the vector tj = (<fii,<f>i) , and using Eqs. JTHJ>, (jT§j) gives 

a fa 

Similarly to what we have found in Sec . 1111 BT1 considering ^-dependence of the one-loop contribution to the effective 
metric, the matrix matter propagator is cancelled upon contraction with the vertex factor. As in the case of the one- 
loop post-Newtonian contributions, the h° terms in the correlation function are associated with the virtual matter 
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quanta propagating near their mass shells. Hence, the first term in the gauge variation of the propagator, Eq. I|74|l . 
results to terms of higher order in the Planck constant. As to the second term, it gives rise to a H° order variation of 
the correlation function, but its structure with respect to the indices \iv is that of an ordinary gauge transformation, 
so 



with some infinitesimal f2 7iQ ^, which proves gauge independence of gauge-invariant functionals built from C^ Q a. Note 
also that as far as variations of the weighting matrix £ a/3 are considered, not only these functionals, but also G^ a/3 
itself turns out to be gauge independent. Indeed, a variation 5£ a P of the weighting matrix is equivalent to the variation 
of the gauge condition F a : 

SF a = 9%Fp , 

provided that 

Then using Eq. J73J) in Eq. (J72J gives 

which implies that C^°J Q/3 is independent of the choice of the matrix C, a f3 ■ 



B. Evaluation of the leading contribution 



Let us proceed to the calculation of the leading contribution to the correlation function, assuming as in Sec. IIII B 21 
that the field producing particle is a scalar described by the action l|45fl . This contribution is contained in the sum of 
diagrams (b)-(f ) in Fig. which has the symbolic form 

where the superscript "tr" means transposition of the indices and spacetime arguments referring to the points of 
observation: [iv «-> a(3, + <-> — , x <-> x' [the transposed contribution is represented by the diagrams collected in part 
(f) of Fig. ©. 

As it follows from the considerations of Sec. II V ATI cf^a can be expressed through C^ a /3 as 



C (0) - lim 



(VT) 2 



d 4 xd 4 x' C^afiix, x') 



{V.T) 



In the DeWitt gauge 



Inuafi reads 



d A zd A z' 



+ D° ++ (x,z) 



D° ++ (x,z) 



D°,(x,z) 



+ D° + _(x,z) 



(z) Vfiv D ++ (z,z') V' a/ 3 <p(z') 

(z)V^ D+-(z,z')V' a/ 3 Uz') 
(z) V> D^+{z,z') V' a0 4>{z') 

(z) V> D—(z,z') V' a p 4>(z') 



D° + _(z',x') 

D°__{z',x') 
D\_{z',x') 

D°__(z',x')\, 



(75) 



(76) 
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where 



Contribution of the third term in the right hand side of Eq. (|76H is zero identically. Indeed, using Eq. l]2Up. and 
performing spacetime integrations we see that the three lines coming, say, into z-vertex are all on the mass shell, 
which is inconsistent with the momentum conservation in the vertex. The remaining terms in Eq. (|76|) take the form 

t < >\ ff d3q d3p a (q) a *(<l + P) ipx'f , \ mn\ 
W(*,») = JJ j^j^-j===ev Po =e q+p -e g , (77) 

where 

__ e ik(x-x ){ m 4 r]tiu7 - la!3 - 2m 2 (ri tiu q a q l3 + T} a ffq^q u ) - Arr^-q^q^kp) 

- 2m 2 T] af) {v(p,(iv) + P(^K) + q^K)) + 4p (A1 <7„) 9(^/3) 

+ Sq^quqakp) + M^k^q^k^ + Aq^qaqp} 
x [D° ++ (k)D ++ (q + k)D° + _(k~p) 

- D° ++ (k)D + _(q + k)D a __(k-p) 

+ D° + _(k)D^(q + k)D°__(k-p)}, (78) 
(M1M2 ' ' ' Mn) denoting symmetrization over indices enclosed in the parentheses, 



• • • Mn) = —. V 

T7 I ^— ' 



Mil M«2 ' ' ' A*in • 

f'1'2 '■'»)- 
porm{12-n} 

As in Sec. IIII B 21 Eq. (|77|l simplifies in the long-range limit to 

= i// ( ^)3 (2^3 IKg)! 2 ^^'"^ W(P, g) , Po = 0. (79) 
The leading term in I^(p, q) has the form [Cf. Eq. (|71Jl ] 

V — P 

This singular at p — > contribution comes from integration over small k in Eq. (|78|l . Therefore, to the leading order, 
the momenta k,p in the vertex factors can be neglected in comparison with q. Thus, 



WM = -^Qa/3 / D° ++ (k)D ++ (q + k)D Q + _(k-p) 



where 



j£k 

- D ++ (k)D + _(q + k)D°__(k- P ) 
+ D Q + _(k)D—(q + k)D Q __(k-p)y (81) 



Furthermore, it is convenient to combine various terms in this expression with the corresponding terms in the trans- 
posed contribution. Noting that the right hand side of Eq. I|81f) is explicitly symmetric with respect to {fJtv} <-> {a/3} 
and that the variables x,x' in the exponent can be freely interchanged because they appear symmetrically in Eq. 1)75(1 . 
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we may write 



(2tt)4 

x { [D a ++ (k)D ++ (q + k)D° + _(k -p) + D°__(k)D^(q + k)D°_ + (k - p)] 
+ [D° + _{k)D—(q + k)D°__(k -p) + D°_ + (k)D ++ (q + k)D° ++ (k - p)] 
- [D ++ (k)D + -( q + k)D a __(k-p) + D a __(k)D- + ( q + k)D a ++ (k-p)]}. (82) 

With the help of the relation 

£>__+£>++=£>+_+£>_+, (83) 

which is a consequence of the identity 

T$(x)4>{y) + T4>{x)4>{y) - 4>{x)^{y) - ${y)$(x) = , 

the first term in the integrand can be transformed as 

D° ++ (k)D ++ (q + k)D° + _(k-p) + D „_(k)D__{ q + k)D°_ + {k - p) 
= D° ++ (k)D ++ {q + k)D°(k - p) + D°(k)D__(q + k)D°_ + (k - p) , 

where 

D(k) = L>+-(fc) + £>-+(fc) = 2m8(k 2 - m 2 ) , D Q (k) = £>(fc)L=o . 

Here we used the already mentioned fact that D+_(k)D + q) = for q on the mass shell. Analogously, the 

second term becomes 

D°(k)D—(q + k)D°__{k -p)- D°_ + (k)D—(q + k)D"(k^p) . 

Changing the integration variables k^k+p, q-^q — p, and then p — > —p, noting that the leading term is even 
in the momentum transfer [see Eq. ijSUJ], and that p(x — x') in the exponent can be omitted in the coincidence limit 
(x — > a;'), the sum of the two terms takes the form 

D°(k) [D ++ (q + k)D° ++ (k -p) + + k)D°_ _(k - p)] . 

Similar transformations of the rest of the integrand yield 

D° ++ (k)D + .(q + k)D Q __(k-p) + D Q __{k)D-+(q + k)D° ++ (k - p) 

-> -i [D° ++ (k)D(q + k)D° ++ (k -p) + D°__(k)D(q + k)D°__(k - p)] . (84) 



Substituting these expressions into Eq. (|82|) and using Eq. I|20|l gives 

d 4 k 
(2tt) 3 ' 

x R J 5[(q + k) 2 -m>] ] 

\ [(k - p) 2 + iO] [{q + k) 2 - m 2 + iO] [k 2 + iO] [{k - p) 2 + iO] J ' { ' 

As was discussed in Sec. II V A II the exponent in the integrand in Eq. (|85|l plays the role of an ultraviolet cutoff, 
ensuring convergence of the integral at large fc's. On the other hand, the leading contribution l|80D is determined by 
integrating over k ~ p where it is safe to take the limit x — > x' . Since (x — x') is eventually set equal to zero, one 
can further simplify the k integral by using the dimensional regulator instead of the oscillating exponent. Namely, 
introducing the dimensional regularization of the k integral, one may set x — x' afterwards to obtain 



25(k 2 ) 



(2tt)3 \ [(k - p) 2 + iO] [(q + k) 2 - m 2 + iO] 
[k 2 + iO][(k - p) 2 + iO] 
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where /i is an arbitrary mass parameter, and e = 4 — d, d being the dimensionality of spacetime. 
Next, going over to the a-representation, the first term in the integrand may be parameterized as 



8{k 2 



5{k 2 



2iri 2 

Substituting this into Eq. 1)86(1 . and using the formulas 



[(k - p) 2 + iO] [{q + k) 2 - m 2 + iO] (p 2 - 2kp + i0){2kq + iO) 
111 dxdydz (e lxk2 + e~" fc2 ) e ^(P 2 -2fcp+*» e42 (2 fcg +*» 



d^ke^ ak +m ^sign(a)- 



f / 7T 



i \\a\ 



2-e/2 







exp 


{- 




\ia 


7re\ 


H 


2i) 


a 



one finds 



K( P ) = fJL 



d 4 ~ e k 



5{k 2 



(2tt) 3 [(k ~p) 2 + iO] [(q + k) 2 - m 2 
ifi e a J^ 3 . . . /f\2-e/2 



(27T) 



■ COS 



dxdydz (^— J 

(?) 



exp 



r l- 



dj/dz e 



-{yp- zqf 



2 z<?) 2 r /2 



exp 



-{yp- zqf 



8tt 2 ' v 4 y * r 27 yy ( OT - zq ) 

Changing the integration variable z — > yz, and taking into account that q 2 = m 2 , qp — —p 2 /2 yields 

\z 2 - (1 + z)a\ t/2 



(pm) e ir c ' 2 
K(p) = — — - — - — cos 



(7re i7r ) 



o 

m 

-e/2 



(?) 



r(l-- / dye*" y*- 1 / dz> , . 



|2 



cos 



(?) r o-D^r- 



|z 2 - (f + z)a 



e/2 



87r 2 m 2 \-p 2 J ' '" \ 4 J " V 2J~ K "'J "' z 2 -(l + z)a 
where a = —p 2 /m 2 . Similar manipulations with the second term in Eq. 1(86(1 give 



L(p) = fx 



d 4 - e k 5[{q + k) 2 -m 2 



dxdydz 



(2nf [k 2 + iO][(k - p) 2 + iO] (2tt) 4 777 

d 4_C fc ^ e "( fc2 + 2fe 9) + e -«(fe 2 +2fcg)^ e iy(p 2 -2fc(g+p)+ l 0) e «(-2fcg+i0) 

^ e rrr°° . . . nx\ 2 -^ 2 



(IX 

dxdydz ( — 



e w <;exp 



--((x-z)g-y(g + p))' 



exp 



-((x + z)q + y(q + p)) 2 



-e/2 



16n 2 m 2 (m) ( 2 

-i7re/4 



'/q ^ U(2/ + ^ + l) 2 +yH 1+e/2 ' [(!/ + z-l) 2 +H 1+e/2 



Changing the integration variables y = ut, z = (1 — u)t yields 

^/ 2 



167r 2 m 2 



(£) r ( 1+ I) 



(87) 



/o Jo \[(t+\) 2 +u(\-u)t 2 a] l +t/ 2 ' [{t-\) 2 + u(l-u)t 2 aY+t/ 2 

On the other hand, L(p) must be real because the poles of the functions D ++ (k), D ++ (k — p) actually do not 
contribute. Hence, 



e/2 



87r 2 m 2 



(£)'~(T) r K) 



dt 



t du 



o Jo [{t+\) 2 +u{\-u)t 2 aY+t/ 2 



24 



Let us turn to investigation of singularities of the expressions obtained when e — > 0. Evidently, both K(p) and L{p) 
contain single poles 



k < \ = 1 r dz 1 1/a i v /T + I 7" + 1 

Adlv W SttW, J Z 2 _ (1 + z)a 87r 2 TO 2 e ^/YTI/^ - 1 ' 1 j 

r M 1 f 1 du 1 1/a Vl + 4/a+l 

dlvW 87r 2 m 2 e 7 [l + «(l-ti)a] 47T 2 m 2 e yiTI^ * ^1 + 4^ - 1 ' 1 ) 

Note that Ldiv = — 2ifdiv Upon substitution into Eq. 1)86(1 the pole terms cancel each other. Thus, C^ va p(p,q) turns 
out to be finite in the limit e — > 0. Taking into account also that divergences of the remaining two diagrams (g), (h) 
in Fig. are independent of the momentum transfer, 10 we conclude that the singular part of the correlation function 
is completely local. 

It is worth of mentioning that although the quantity Kdi v {p) is non-polynomial with respect to p 2 , signifying non- 
locality of the corresponding contribution to Cnua/3( x > x ')i it is analytic at p = 0, which implies that its Fourier 
transform is local to any finite order of the long-range expansion. Indeed, expansion of K^ivip) around a — reads 



a 



(-l) n (nl) 2 



lD7r z m z e \ 6 (2n + lj! / 

Contribution of such term to C^ v {x,x') is proportional to 

/ w?^*"-*" ( l - S + - •) = sm <*' - - J + ^ ,3, <*' -*)+■••■ 

The delta function arose here because we neglected spacial spreading of the wave packet. Otherwise, we would have 
obtained 



d 3 q d 3 p 

(2tt) 3 (2tt) 3 ' 



&*(g)% + p)e- 4p(x '- xo) (l-^- + ---) =0 for cc' -Ko^lU: 



as a consequence of the condition 1)50(1 . In particular, applying this to the correlation function, we see that its divergent 
part does not contribute outside of W. Thus, the two-point correlation function has a well defined coincidence limit 
everywhere except the region of particle localization. 

Perhaps, it is worth to stress once more that the issue of locality of divergences is only a technical aspect of our 
considerations. This locality does make the structure of the long-range expansion transparent and comparatively 
simple. However, even if the divergence were nonlocal this would not present a principal difficulty. A physically 
sensible definition of an observable quantity always includes averaging over a finite spacetime domain, while the 
singularity of the two-point function, occurring in the coincidence limit, is integrable (see Sec. IIV A 1(1 . The only 
problem with the nonlocal divergence would be impossibility to take the limit of vanishing size of the spacetime 
domain [Cf. Eq. J75J)]. 

Turning to the calculation of the finite part of K(p), we subtract the divergence l(89|) from the right hand side of 
Eq. I(88|) . and set e = afterwards: 



K Rn (p) = lim[K(p)-K div (p)] 



1 



16ir 2 m 2 



Tn 7r + 7 + 2 In 



fj,m 



dz 



z 2 - (1 + z)a 



00 lnlz 2 - (l + z)a\ I , . 

o z 2 -(l + z)a 



10 The corresponding integrals do not involve dimensional parameters other than p, and therefore are functions of p 2 only. Hence, on 
dimensional grounds, both diagrams are proportional to 

where ci 2 are some finite constants. 
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where 7 is the Euler constant. The leading contribution is contained in the last integral. Extracting it with the help 
of Eq. (|A6() of the Appendix, we find 

#W(p) = _J— . (93) 
o4my — p 

As to L(p), it does not contain the root singularity. Indeed, 

r~ e/2 ( f ne \ r ( 1 + e / 2 ) 



V ' 8ir 2 m 2 Km) K 4 



and therefore £fi n (p) = lim[L(p) — ^div(p)] is finite at p = 0. It is not difficult to verify that L^ n (p) is in fact analytic 

e— >0 

at p = 0. Thus, substituting Eq. into Eqs. JSHJ, (|7§1) . and then into Eq. {7SJ, using Eq. (J5SJ, and restoring 
ordinary units, we finally arrive at the following expression for the leading long-range contribution to the correlation 
function 

e£U*) = S ^ 2J ^ > r = |x - xo| . (94) 

This result coincides with that obtained by the author in the framework of the S-matrix approach |27| . The nonrel- 
ativistic gravitational potential <& s is related to the 00-component of metric as $ s = hooc 2 /2 . Hence, the root mean 
square fluctuation of the Newton potential turns out to be 



(A*<> = f§. (95) 

Note also that the relative value of the fluctuation is 1/ Y2. It is interesting to compare this value with that obtained 
for vacuum fluctuations. As was shown in Ref. [2]]], the latter is equal to y/2 (this is the square root of the relative 
variance A 2 used in Ref. HJ). 

We can now ask for conditions to be imposed on a system in order to justify neglecting quantum fluctuations of 
its gravitational field. Such a condition can easily be found out by examining dependence of the h° contribution on 
the number of field producing particles. The structure of <^-lines in the diagram of Fig. |H{b) is the same as that 
in Fig. Ela). Therefore, in the case when the gravitational field is produced by a iV-particle body, this diagram is 
proportional to m 2 N, where m is the mass of a constituent particle. Correspondingly, the root mean square fluctuation 
of the potential is proportional to m\/~N — M/y/N (M is the total mass of the body), while its relative value, to 
l/VN . We thus arrive to the condition of macroscopicity of the system, found in Sec. IIII B~2l in connection with the 
post-Newtonian loop contributions. In the present context, it coincides with the well-known criterion for neglecting 
quantum gravity fluctuations, formulated in a somewhat different manner in Refs. pH l28l l2fll| . 

C. Orbit precession in the field of black hole 

Here the results obtained in the preceding sections will be applied to the investigation of dynamics of a classical 
particle in the gravitational field of a black hole with mass M. The particle will be taken testing, i.e., it will be 
assumed sufficiently light to neglect its contribution to the gravitational field, and sufficiently small compared with 
r g = 2GM/c 2 to allow considering it as a pointlike object. The latter assumption, in particular, justifies the use of 
the expression l|94|l for the correlation function, obtained in the coincidence limit. 

In order to find equations of motion of the particle, we have to calculate its effective action. The action functional 
for a pointlike particle has the form 



S t = -m t J V 1 g^ v dx^dx v = -m t J dr^/l 



— i ■ a dxtl 

+ hu V x^x v , a; M = 

dr 



where dr is the particle proper time defined with respect to the flat metric, dr 2 = n^dx^dx 1 ' . Inserting this expression 
into the generating functional integral of Green functions, one can write, in view of the assumed smallness of the test 
particle contribution, 



J 2)$_y" D$ +ex p{i(5 + S r t +i^)} 

= ZskIa^o + Jvt>_ J V^+iSt exp{i[S + j<t>)} 



?SK(l + *<in|5 t |in))|. Aj=0 . (96) 

K = 
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The generating functional of connected Green functions takes the form 

W = -tln{Z SK (1 + i(in|S t |in))} = W SK + <in|5 t |in) . 
Taking its Legendre transform, we see that the effective particle action is just 

T t = (mlStlin). 

Expanding St in powers of h llVl the right hand side of this equation can be evaluated in the first post-Newtonian 
approximation as 



-ra t 



dr\l + \{K v )x»x v - l - [<V)M + tfLp 



x^x v x a x fi 



-rot J \j (v^ + (h^))dx»dx" ~ ie^^ ) 2 . 



(97) 



Substituting Eqs. (|(j6fl . (|94|l for the mean gravitational field and its fluctuation, we see from Eq. 1|97[1 that the test 
particle motion in the fluctuating field of a black hole is effectively the same as in a non-fluctuating gravitational field 
described by the following spacetime metric 

ds 2 = gf v dx»dx v = \ 1 - T JL - i ) dt 2 - ( 1 + H - ] dr 2 



-r 2 1 



7ri 



JL _i_ g 



8r 2 



(de 2 + sin 2 9 dtp 2 ), 



4r 



(98) 



Now, it is not difficult to calculate the secular precession of the particle's orbit. Let the test particle with move 
in the equatorial plane (6 = it/2) around black hole. Then we have the following Hamilton- Jacobi equation for the 
action T t 



llv dT t 8T t 



9cs 



dx^ dx u 

where is the reciprocal of <?°^. A simple calculation gives, to the leading order, 
S h = -Et + Lip 



J dr [(£' 2 + 2mtE') + ^ (m 2 + 4m t E>) - 1 (l 2 - ^r 2 m t 2 ) 



1/2 



(99) 



where E, L are the energy and angular momentum of the particle, respectively, and E' = E — m t its non-relativistic 
energy. The first two terms in the integrand in Eq. (j99(l coincide with the corresponding terms of classical theory, 
while the third does not, leading to the angular shift of the orbit 



Sip = 



UirGM 



2c 2 a(l-e 2 ) 

per period (a and e are the major semiaxis and the eccentricity of the orbit, respectively). 
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V. CONCLUDING REMARKS 



In this Chapter, we have determined the long-range behavior of the gravitational fields produced by quantized 
matter in the first post-Newtonian approximation. From the point of view of establishing the correct correspondence 
between classical and quantum theories of gravitation, Eq. (|56|l demonstrating violation of general covariance by the 
fi°-order loop contributions, and Eq. l|MJl describing ^-fluctuations of the spacetime metric turned out to be of special 
importance, and led us independently to a macroscopic formulation of the correspondence principle. In turn, this 
formulation endows H° loop contributions with direct physical meaning as describing deviations of the spacetime metric 
from classical solutions of the Einstein equations in the case of finite number of elementary particles constituting the 
gravitating body. 
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In connection with these results, it is worth also to make the following general comments. 

As we have seen, an essentially quantum character of elementary particle interactions makes classical consideration 
inapplicable to systems whose dynamics is governed by interactions of a finite number of constituent particles. On 
the other hand, there is a deep-rooted belief in the literature that the quantum field description of interacting remote 
systems, each of which consists of many particles, is equivalent to that in which these systems are replaced by 
elementary particles with masses and charges equal to the total masses and charges of the systems. In other words, 
without calling it into question, the familiar notion of a point particle is carried over from classical mechanics to 
quantum field theory. This point of view is adhered, for instance, in the classic paper by Iwasaki [2j where it is 
applied to the solar system to calculate the secular precession of Mercury's orbit, considered as a "Lamb shift". The 
Sun and Mercury are regarded in Ref. as scalar particles. As we saw in Sec. II V Bl under this assumption the 
root mean square fluctuation of the solar gravitational potential would be 71% of its mean value. Fortunately, such 
fluctuations are not observed in reality. This is because the Sun is composed of a huge number of elementary particles 
each of which contributes to the total gravitational field. As we have seen, the relative quantum fluctuation turns out 
to be suppressed in this case by the factor 1/y/N ~ yj m P roton/-Mo ~ 10~ 28 . 

Another example of attempts to recover the nonlinearity of a classical theory through the radiative corrections can 
be found in Ref. [3(| • The authors of [3(j claim that the electromagnetic corrections of the order e 2 to the classical 
Reissner-Nordstrom solution are reproduced by the diagram of Fig. |5Ja) in which the internal wavy lines correspond 
to the virtual photons. However, as we have shown, it is meaningless to try to establish the correspondence between 
classical and quantum theories in terms of elementary particles, because quantum fluctuations of the electromagnetic 
and gravitational fields produced by such particles are of the order of the fields themselves. On the other hand, 
because of its inappropriate dependence on the number of particles, the diagram |2fa) fails to reproduce the classical 
physics in the macroscopic limit. This can be shown using the same argument as in the case of purely gravitational 
interaction. Namely, given a body with the total electric charge Q, consisting of N = Q/q particles with charge q, the 
contribution of the diagram |2Ia) is proportional to N ■ q 2 — Q 2 /N turning into zero in the macroscopic limit. The 
relevant contribution correctly reproducing the e 2 -correction to the Reissner-Nordstrom solution is given instead by 
the tree diagrams like that in Fig. ^b) in which internal wavy lines correspond to virtual photons. 

Finally, we mention that investigation of quantum fluctuations of the Coulomb potential similar to that carried out 
in Sec. IIVI shows that the macroscopic formulation of the correspondence principle given above extends also to the 
case of electromagnetic interaction [22J |31| • 



where A, B > are some numbers eventually set equal to 1. Performing an elementary integration over u, we get 



APPENDIX A: ROOT SINGULARITIES OF FEYNMAN INTEGRALS 



The integrals 




encountered in Sec. IIII B 21 can be evaluated as follows. Consider the auxiliary quantity 





Now consider the integral 




(Al) 



taken over the contour C shown in Fig. [7| J(A, B) is zero identically. On the other hand, 
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W± denoting the poles of the function f(w, A, B). 



A B A 2 

Change z — > — A — z in the first integral. A simple calculation then gives 

J(A,B) = - 7 =B-V 2 [1 + —,) - dt — * > ^ (A2) 



2^a V 4B / 2 Jo B + at(A-t) 

The roots are contained entirely in the first term on the right of Eq. I|A2(I . The integrals J nm are found by repeated 
differentiation of Eq. I|A2J| with respect to A, B. Expanding (1 + aA 2 / AB)^ 1 / 2 in powers of a, we find the leading 
terms needed in Sec. IIII B 21 

7T 2 StT 2 

J™ ot = - — VS, J™ ot = -^!^, root^^lSvr 2 ^ (A3) 
31 16 32 32 33 128 

Next, in the course of evaluation of the integral K(p) we encountered the integral 

ln^-(l+goj = r rfz ln|( Z - Pl )(,- P2 )| 



z 2 -(l + z)a y o («-pi)(«-^2) 

a ± \J a 2 + 4a 

Pl,2 = g • 

To evaluate this integral, let us consider an auxiliary integral 

A — \ dw M(™-pi)(™-w)] (A4) 

taken over the contour C in the complex w plane, shown in Fig. [S] A is zero identically. On the other hand, for 
sufficiently small positive a, one has 

. f pi - a , \n[{z - pi){z - p 2 )] . In(p2-pi) i7rlna + 7r 2 /2 

A = I UZ ; 77 : Z7T- 



(z-p 1 )(z-p 2 ) Pi- Pi Pi- Pi 

P2 ~ a ln[(z — pi)(p2 — z)] — in . \n(p 2 — Pi) — in in In a + 7r 2 /2 
(z-pi)(z-p 2 ) Pi- Pi Pi- Pi 



pi+a 

ln[(z -pi)(z- pi)] - 2in 



lp 2 +a " (z-Pl)(z~Pl) 

Changing the integration variable z — > (pi + p 2 ) — z in the first integral, and rearranging yields in the limit <r — > 

a = i f P1+P2 dz H^-^-p;)] + je!ZL . (A5) 



The root singularity is contained in the second term, because 

(z-pi)(z-p 2 ) 



ln[(p 2 -z)(z-p 1 )] 
rfz — - 77 ; ► — In a for a — ► 0. 



Thus, 



^ 2 



^root = (A6) 
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(a) (b) 

FIG. 1: Lower-order tree diagrams representing solution of the Einstein equations, (a) The first order (Newtonian) term, (b) 
One of the second-order diagrams describing the first post-Newtonian correction. Wavy lines represent gravitons, solid lines 
the auxiliary source t^" , or expectation value of the energy-momentum tensor of matter, {T y ' v ). 




FIG. 2: Some of the one-loop diagrams contributing to the effective gravitational field, (a) The only diagram giving rise to the 
root singularity with respect to the momentum transfer p. (b) An example of a diagram free of the root singularity. Wavy lines 
represent gravitons, solid lines massive particles, q is the particle 4-momentum. 



+ + - - + 




FIG. 3: Some of diagrams arising upon writing out the matrix diagram in Fig.|3a) according to the Schwinger-Keldysh rules, 
(a) The only diagram giving rise to a nonzero contribution in the long-range limit, (b) Diagram vanishing in the long-range 
limit because of the condition p° = 0. (c) An example of a diagram which is zero identically because of the vanishing of one of 
its vertices (the left cf> 2 h vertex). 




FIG. 5: Diagrams responsible for the nontrivial contribution to the right hand side of Eq. I|47|l . 
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FIG. 7: Contour of integration in Eq. HA It . 
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FIG. 8: Contour of integration in Eq. 1A4I . 



